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ABSTRACT 

Flaring release of magnetic energy in solar corona is only possible if the magnetic 
field deviates from a potential one. We show that the linear MHD modes excited 
on top of the non-potential magnetic field possess a nonzero kinetic helicity. Ac- 
cordingly, this necessarily results in a noticeable kinetic helicity of the turbulence, 
composed of these linear modes with various scales and random phases, gener- 
ated at the flare site by the primary energy release, which may be important for 
many applications. In particular, a nonzero turbulence helicity has a potentially 
strong effect on the particle acceleration because the helical component of the 
turbulence induces a mean regular large-scale (DC) electric field capable of di- 
rectly accelerating the charged particles in addition to the commonly considered 
stochastic turbulent electric field. In this paper, we derive the kinetic helicity 
density of the linear MHD modes excited on top of a twisted large-scale mag- 
netic field, estimate the corresponding turbulence helicity and take its effect on 
stochastic particle acceleration by the turbulence into consideration; in particular, 
we compare this induced mean electric field with the electron and estimated ef- 
fective ion Dreicer fields. We have discovered that this, so far missing but highly 
important, ingredient of the turbulence at the flare site can be responsible for 
the thermal-to-nonthermal energy partition in flares by controlling the process of 
particle extraction from the thermal pool and formation of the seed particle pop- 
ulation to be then stochastically accelerated to higher energies. In addition, it is 
naturally consistent with such puzzling flare manifestations as spatial separation 
of electron and proton emission sites, electron beam formation, and enrichment 
of the accelerated particle population by 3 He and other rare ions. 

Key words: Sun: flares — acceleration of particles — turbulence — diffusion — magnetic 
fields — Sun: radio radiation 



1 INTRODUCTION 



Stochastic acceleration by turbulence is currently believed to pla y a dominant role in charged particle acceleration 
in impulsive solar flares. Ind eed, it has been recognized long ago ( Hamilton fc Petrosian 19921 : Petrosian et al. 1994 : 
Miller et all 19971 : lMillerlll99cf ) that models of stochastic acceleration of charged particles by a postulated MHD or whistler 
wave ensemble are overall consistent with global manifestations of the acceleration process such as total energy content 
of the accelerated particles, the spectrum shape, the highest energy of the accelerated particles, and the acceleration 
time. Furt hermore, the stochastic acceleration was shown to be consistent with a number of more detailed acceleration 
properties ijCrigis fc Benzlbood : uBykov fc Fleishmanll2009h e.g.. soft-hard-soft (SHS) spectrum evolution established in 
detail based on t he RHESSI observatio ns jGrigis fc Benzn2004r ). In addition, radio detection of the very acceleration 
region in a flare (jFleishman et alj|201ll ) is suggestive in favor of a stochastic acceleration process with roughly energy 



© 2012 RAS 



2 Gregory D. Fleishman and Igor N. Toptygin 



independent el ectron residence time in the accelerat ion region, see also recent review papers on the particle acceleration 
in solar flares (|Zharkova et alj|201ll ; |Petrosianll2012r ). 

Nevertheless, the standard paradigm of the stochastic acceleration is not apparently the whole story since it cannot 
easily predict the acceleration efficiency leading to the observed variety of the therma l-to-nonthermal energ y partitions 
from 'heating without acceleration' (at least, a t the early stage of some flares, see e.g.. iBattaglia et al.ll2009f ) to 'acceler- 
ation without heating' (jFleishman et al.ll2011T ). which, perhaps, requires a specific mechanism capable of dissipating the 
flare energy into the plasma heating and at the same tim e extracting a fraction of particles from the thermal pool and 
supplying them to the stochastic acceleration process, see I Altvntsev et al.l (|2012r i for more detail. Furthermore, it is at 
odds , e.g., with particle accelera tion in form of beam^J (detected via radio type III bursts or microwave continuum, see, 
e.g., Altvntsev et al. 2007 . 2008h . It is no t easily consiste nt with huge enhancement of some ion abundances including 
3 He and ultra heavy ions (see, however, Liu et al-lboOrj) observed in solar energe tic particle (SEP) events produced 



from impulsive solar flares (|Kocharov fc Kocharqv 



enhanced positron production at the flare site (|Kozlovskv et al 



1984; iReames 



1993; lMasonll2007l ), and also indirectly impl ied by the 



20041 ). In addition, [jhare fc Murphvl |l99St ) analyzed 



various gamma-ray lines produced by different accelerated ions at the flare site and found that the He/p ratio was 
noticeably enhanced (by a factor of 5) relative to the standard photospheric value. 

Then, the RHESSI discovery of distinctly different spatial lo cations for hard X-ray emission produced by accelerated 
electrons and gamma emission produced by ions in the same flare ( Lin et ai]|2003l ; Hurford et al.|[200r3 ; IVilmer et alj|2011 ) 
puts further constraints on particle acceleration in flares. The first of the events, 23 July 2003, where this separation was 
observed, displayed a complex source structure with the HXR (gamma-) source associated with shorter (lar ger) post-flare 



loops , which has been interpreted in terms of resonant stochastic acceleration by cascading MHD turbulence IjEmslie et al 



2004). Since then, a few more strong gamma-ray flar es have been studie d. Although no statistically significant set of 



flares with gamma-ray imaging has been accumulated (jVilmer et al.ll201lr ). the several available events display a variety 
of spatial relationships with a significant spatial separation between the electron and ion emission sites. In t he cases (tw o 



of five available events) where the ion emission site coincided with a HXR footpoint (which seems at odds to lEmslie et al 



(2004) model), it was a weaker HXR footpoint, indicating that the accelerated electrons and ions preferably precipitate 
into the opposite footpoints of a flaring magnetic loop or loop system; something that could also happen in the other 
events, although the magnetic connectivity between sources is less clear in those cases. 

These observed properties of particle acceleration in flares imply that an important ingredient is miss ing from our 



understanding of the flare acceleration process commonly ascribed to the stochas tic acceleration mechanism (jMiller et al 



19971 ; lAschwandenl2002l ; lEmslie et ai1l2004l ; Izharkova etahlboill ; |Petrosianll2012l ) . Indeed, within the concept of stochastic 



acceleration, in contrast to direct acceleration by a DC electric field, the particles being scattered by stochastically 
turbulent magnetic fields move diffusively (at least, inside the acceleration region) and, according to the Fisk law, form 
a particle flux in the direction opposite to their density gradient, which is supposed to be independent on the particle 
charge sign. Some asymmetry between two opposite footpoint emission is expected due to asymmetry of the magnetic 
loop; however, this asymmetry must be the same for both protons and electrons, not the opposite as observed. Having 
the protons and electrons precipitating into the opposite footpoints would be natural in presence of a DC electric field. 
The DC electric fields can accelerate the charged par ticles directly but to a modest energy ( ^ 100 keV) only (e.g., 
Holman|[l985l ; lAschwand cn 2002; IZharkova et alj|201ll . and references therein); thus, an acceleration process combining 



features of both stochastic acceleration and acceleration in DC electric field is called for. 

We are not going, however, to arbitrarily postulate these two highly distinct ingredients — a regular DC electric field 
and random turbulent fields needed for stochastic acceleration. Instead, we envision a model where the appearance of 
regular electric fields is an intrinsic property of the very same turbulence that gives rise to the stochastic acceleration — 
a helical turbulence model — in which the effect of helicity does not substitute but complement the more well-known 
stochastic acceleration process. We demonstrate below that the helicity of the turbulence created in flaring loops is an 
unavoidable direct consequence of the magnetic field nonpotentiality at the flare site, which is itself needed to drive 
the flare energy release. Furthermore, we show that the turbulence kinetic helicity in the corona is coherently built up 
from non-zero helicity of individual linear MHD waves forming the turbulence, which is linked to the current helicity 
and so can reliably be estimated using measurements of the current helicity at the photosphere. Stated another way, 
the helicity is an unavoidable property of the turbulence created in the nonpotential, twisted flaring loops; account of 
the turbulence helicity effect on the stochastic acceleration offers, as we show below, a nice natural way of addressing 
the mentioned outstanding questions about particle acceleration in flares including the thermal-to-nonthermal energy 
partition, the electron/ion spatial separation, formation of electron beams, and significant 3 He enrichment. 



1 Alt hough the beams can well be formed outside the acceleration region due to transport effects fe.g,. lReid et al]|201ll ; Ivilmerl 
I2012T) 
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2 PARTICLE ACCELERATION BY HELICAL TURBULENCE 



Theory of particle acceleration by helical turbu l ence was propo sed long ago bv lKichatinov (1983) and then developed by 



Fedorov et all (|l992h ; lFedorov fc Stehlikl (|201oh ; iFedorovl (|201ll . and references therein) in relation to Galactic cosmic ray 



acceleration problem. Since it has not been discussed yet in the solar context, we outline below the essence of the theory. 
Let us denote the exact (fluctuating) distribution function of the fast particles of a given kind as f(r,p,t); apparently, 
it satisfies the collisionless kinetic equatioijf]: 

The force 

F= eE+ -v xB =~(v-u) x B (2) 

c c 

contains both electric and magnetic components, where we employ the high plasma conductivity providing E — —uxB/c. 

Let us consider the Lorenz force structure in the presence of a regular magnetic field Bo(r) with a turbulent wave 
spectrum superimposed. The turbulence can be characterized by a random magnetic field vector b and a random vector 
of the fluid velocity u (the fluid is adopted to be at rest so no regular flow velocity is present). Thus, substituting 
B = Bo + b and u into Eqn ((2j) we obtain 

g 

T — - [v x Bo + v x b — u x Bo — u x b] . (3) 

The first term describes the particle gyration in the regular magnetic field Bo, the second one — the particle angular 
scattering by magnetic irregularities b, while the third one — the standard stochastic (second order Fermi) acceleration 
by the turbulent electric field E — —u x Bo/c. Remarkably, the second and third terms have zero means; thus, they 
result in the second order diffusion in the momentum space — isotropization and stochastic acceleration, respectively. The 
fourth term, also representing an electric field — u x b/c is, however, different from the third one because it can have 
a non-zero mean (it x b) 7^ 0, which will represent a regular electric field capable of directly accelerating the charged 
particles. Importantly, having this nonzero regular electric field constrains the turbulence in a certain way. Indeed, based 
on the symmetry consideration, the mean electric field E = (u x b) /c must be directed along the mean magnetic field, 

E = pB , (4) 

and given that the electric field is polar, while the magnetic field is axial vector, we conclude that the coefficient p must 
be a pseudoscalar. Stated another way, only such turbulence whose measures intrinsically include a pseudoscalar can 
result in a nonzero regular electric field cE = (u X b). In what follows we show that this is the case of helical turbulence. 
To consider this effect quantitatively let us write down Eqn. fl} in the form: 

f + + B.B/ = 0, (5) 

where 

D= -(v-u) x A (6) 
c dp 

is the operator that changes the value and direction of the particle momentum. As has been said above we adopt u to 
be a purely random vector with zero mean (u) — 0. 

Let us average Eqn ((SJ) over the turbulence ensemble as is commonly done to address the stochastic acceleration 
problem. For simplicity, we assume that the average distribution function does not depend on coordinates, so that 

(f(r,p,t))=F(p,t), f = F + Sf(r,p,t). (7) 

After averaging we obtain the equation 

where each term is primarily linked with a certain component of the Lorenz force in Eqn ([3}: the first term at the rhs 
originates from averaging (bSf(r,p,t)) and describes angular scattering of the particles by magnetic irregularities, which 
does not change the particle energy, while the second term originating from the electric field averaging, (uSf(r,p,t)), 
describes the standard stochastic acceleration by turbulence. The angular scattering term is expressed via the particle 



2 The collisional integral is discarded assuming the acceleration dominates over the Coulomb losses; we return to the role of 
Coulomb collision, that is important for thermal particles, later in §[4] 
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mfp A(p) that is set up by the turbulent magnetic fluctuation spectrum. To be specific, we adopt the power-law spectrum 

(iz-ijfcr 1 ^) 



(F) k dk 



k" 



-dk, 



v > 1, k c < k < 2n/l n 



(9) 



where 6 = ^/ (6 2 ) is the total rms turbulent magnetic field so that lb 2 ) = J (b 2 )kdk, v is the turbulence power-law 
spectral index, L c = 2ii/k c is the main (largest, energy-containing), while / m i n is the minimal scale of the turbulence, 
which can decrease at the course of the turbulence cascading in contrast to the main scale L c , which remains roughly 
const ant. Then, the mfp A is given by (§ 15.3, pp 207-211, in lToptvgin|[l98i ; ISchlickeiseij 1 19891 : lHamilton fc Petrosianl 
19921 . Eqns 78, 13, respectively) 

2-v 



A 



3. 



oc p 



(10) 



r) 

O = v x T^jj is an angle-alternating operator, and the operator S a includes all acceleration effects other than that 
produced by the nonzero mean electric field E = {u x b) jc if present. 

Below we concentrate on this new effect described by the term proportional to ( !F) , qualitatively distinct from what 
is commonly included in the standard Fermi acceleration. As has been noticed based on qualitative consideration of the 
Lorenz force in form ((3|, if u and b at a given point are correlated and not parallel to each other then {u x b) 7^ 0, which 
yields 

(F) = ~vxB --(uxb), (11) 
c c 

where the last term represents an electric force {J-h) acting on the particles. 

In what follows we employ a d i scovery made by Parke r, Steenbek and Krause (see IVainshtein fc Zerdovichlll972l ; 
Parkerj|l979l : IVainshtein et al.lll980l ; iKrause fc Radlerlll980r ) within the turbulent dynamo theory that large-scale fields 
can be produced by a turbulence lacking mirror symmetry (helical turbulence). This ability to amplify the large-scale 
magnetic field is the main driver of studying the helical turbulence in the astrophysical context, although for the particle 
acceleration problem addressed here we are primarily interested in the accompanying large-scale electric field, which is 
an unavoidable byproduct of the turbulence kinetic helicity; see Appendix II. 

The correlation tensor of this turbulence contains both polar and axial vectors. Perhaps the simplest example of 
such structure represents the correlation tensor of statistically uniform gyrotropic incompressible fluid: 



U^{k,t) = A{k,t) S, 



k 2 



iP(k,t)e^xk\. 



(12) 



Importantly, the function P(k,t) here must be a pseudoscalar (recall, the presence of a pseudoscalar is a necessary 
condition for having a nonzero mean electric field E), i.e., it changes the sign at the coordinate inversion. The turbulence 
helicity can be quantified by the following integral parameter 



1 r°° 

- {u(r,t + r) -V x u(r,t))dT = 
Jo 



T -^{u(r,t) ■ V x u(r,t)} 



3tt 2 



k i P{k,0)dk. 



(13) 



This measure, a, is called the kinetic helicity parameter, which is defined by the second term in Eqn (|12[) only. 

Apparently, a nonzero helicity is only possible for the systems where a pseudoscalar can be built from available 
physical parameters; for example it is the case of the dot product of a polar vector (e.g., gradient of pressure or density, 
which requires plasma inhomogeneity) and an axial vector (e.g., the magnetic field). In particular, in the magnetized 
coronal plasma the helicity can well be driven by the twist of the nonpotential field lines, where the polar vector of the 
electric current density j is available. It is, therefore, important that the mean magnetic field is itself twisted, so that 

B = B + b, (B) =B ^0, VxB,/ 0, (b) = 0, (14) 

where Bo(r) is a mean regular nonpotential magnetic field in the acceleration region and b is the fluctuating turbulent 
magnetic field. 

As soon as the helicity parameter has been specified (we return to estimating it later, see Appendixes), the averaged 



Jia OWV-ZIJ. <jLO ^7 ^-'Llji LUlllV.- L VjI IILUO Ul^Vjll □p^VlllV-U ^ YVIj 1 lu U 111 11 LU VjlJUIIIICHjIII^ J. Vj ±UL\jKjL -j Okjkj J 1 ^ j ljl±X_- <_L V «^ 1 U.^LjVI 

electric force is straightforwardly calculated w ithin the turbulent dynamo theory I Vainshtein fc Zerdovichl 19721 : Parker 
19791 : IKrause fc Radlerill980l : iKichatinovll 19831 ) and neglecting the dissipation (and magnetic field nonuniformity) it has 
the form agreeing with expectation Q 



'■(u x b) — eEh 



-aB , 



(15) 



which depends on the turbulence helicity only, while does not depend on other plasma properties, e.g., on the plasma 
beta. Hence, Eqn Q for the particle distribution function receives the form 
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OF ec „ si „ ea „ dF v ^2 

-^ + J B .OF--B .- = WF) O F, (16) 

where for now we have neglected the contribution from the term S a F, since we are interested in the contribution coming 
from the large-scale electric field (-a/c)Bo; however, we will compare efficiencies of these two acceleration processes 
later. 

The regime of particle acceleration including the acceleration rate in the DC electric field depends on the particle 
transport regime, which is controlled by the particle mfp A(p). At a relatively small scale, I <C A, the particle guiding 
center moves almost rectilinearly along the mean magnetic field and the energy gain of a nonrelativistic particle is 
proportional to t 2 , i.e., the energy gain is fast. Apparently, at this stage the particles can form a beam and quickly gain 
energy of the order o£j eZEhA, where eZ is the charge of the particle. At larger scales, I 2> A, the particle experiences 
many scatterings and so moves diffusively back and forth along the magnetic (and, having the same direction, DC 
electric) field. Accordingly, the particle gains energy from the DC electric field when it moves one direction, while loses 
it when moves the opposite direction; thus, the particle energy changes much slower with time, than at the small scale 
regime. Even though the particles move diffusively, the process of their energy gain from the regular DC electric field 
is different compared with that from stochastic electric fields during standard stochastic acceleration process. Indeed, 
consider the energy change of a particle, which comes back to its original location after it has randomly walked for a 
while. In the case of purely stochastic electric fields this change is also stochastic — it can be positive or negative, being 
positive on average. In contrast, in the case of purely regular DC electric field this energy change is zero equivalently. 
Here, the mean positive energy gain in time is provided by regular unidirectional drift of the particles, driven by this 
DC electric field, superimposed on the particle diffusion. Stated another way, regardless of the particle transport regime 
the particle energy gain in a DC electric field is a regular function of its original and final spatial locations; for the small 
scales (small source size) it also a regular function of time, while for large scales (large source size) it is a random function 
of time; thus, individual particles can either gain or lose their energy depending on their individual paths, while gaining 
the energy with time on average even though the accelerating field is a regular DC electric field. 

Since the solution for the small scale regime, I <C A, is trivial, we consider the other regime, I ^> A, in more detail. 
In a general case kinetic equation (|16[1 describes anisotropic particle distributions. However, if the departure from the 
isotropy is small, which is reasonable to expect within the considered transport regime because of the efficient angular 
scattering of the particles by the very same turbulence described by the rhs of Eqn (|16[) . this equation can be further 
simplified by using the diffusion approximation (we emphasize that this will also be very convenient for the purpose of 
further comparison of this mechanism with the temporarily discarded statistical Fermi acceleration) if the source size L 
is much larger than the particle mfp L ^> A. Assuming this to be the casjf], we adopt a weakly anisotropic distribution 
function to have the form: 



F(r, P ,t) = ± 



N{r,p,t) + ^v-J(r,p,t) 

v 



(17) 



where 

N{r,p,t) = j F{r,p,t)dQ p , 

J(r,p,t) = J vF(r,p,t)dQ p (18) 

are, respectively, the isotropic part of the distribution function and the flux density of the particles with a given energy; 
weak anisotropy implies J<^vN. 

Substituting Eqn (|17|) into Eqn (|16|l , we obtain the set of coupled equations 

dN ea h 

c 2 J ^ dp 



-Bo 

at cpv 



(19a) 



1 dJ ec „ , ea dN „ 1 r . 
-- + -B xJ--— B = --J 19b 

v at vt 3c dp A 

for the isotropic part of the particle distribution function N and the flux density J. Note that the direction of the 
flux density depends on the charge sign, implying opposite flux directions for the protons and electrons, respectively. 



3 This estimate is valid if the source size L is larger than the particle mfp; otherwise the estimate apparently reads eZE h L. 

4 For example, adopting a modest turbulence level of b 2 /B 2 ~ 10™ 4 , L c ~ 2 ■ 10 7 cm, v = 1.5, and Ro ~ 0.2 cm, Eqn 110t 
yields A ~ 2 ■ 10 7 cm for electrons, which is indeed much smaller than the size of acceleration region, ~ 10 9 cm, determined from 
observations (e.g.. iFleishman et al.|[201ll V 
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Remarkably, it also depends on the a sign. For a flaring loop in the solar corona this means that for a given magnetic 
polarity of the loop footpoints at the photosphere level, the electric field direction (and, accordingly, directions of positive 
and negative charge flows) differs, respectively, for the left and right twists of the magnetic flux tube — a prediction, which 
can be tested observationally. 

Let us eliminate J from (|19a|l neglecting the term (l/v)8J /dt in Eqn (|19b|) . We then arrive at the diffusi on equation 
for th e isotropic part N of the weakly anisotropic distribution function F, Eqn (|17[) , in the momentum space (jKichatinov 
19831 ), 

— = — — p 2 D h — (20) 
dt p 2 dp dp ' 

2 a 2 2 — v 

a Ap p 

Dh = o..T>2 — ( 21 ) 



where 



3vK 



is the diffusion coefficient which describes the particle acceleration by the helical turbulence, and Ro = cp/eBo is the 
Larmor radius. Eqn (|20[1 ignores finite size of the acceleration region and so permits acceleration up to arbitrary high 
energy; in fact, for a flux tube with a length L the maximum energy can apparently be estimated as £ ma x ~ SeE^L. 

In spite of the presence of the regular electric field along the regular magnetic field, we have arrived at a diffusive 
form (in the momentum space) of equation describing particle acceleration by the helical turbulence. As has been noticed 
above, this happens because of efficient particle angular scattering by the turbulence, which makes local particle motion 
back and forth relative to the accelerating electric field almost equally probable, although the probability to gain energy 
remains larger than to lose energy. Accordingly, though the regular electric field enhances explicitly only the parallel 
particle momentum, the angular scattering transfers it to the transverse momentum, so all momentum components rise 
proportionally roughly preserving the distribution isotropy. 

In fact, the (discarded earlier) second order Fermi-acceleration process does contribute to the right-hand side of Eqn 
( 20l in addition to the acceleration due to turbulence helicity. The Fermi-acceleration diffusion coefficient is given by 



Tbptvgin|[l985l : iFedorov et alll 19921 ) 



Df = oc C (22) 

3vA v y ' 

Thus the helical part of the velocity correlation tensor gives a dominant contribution to the acceleration if 

D h a 2 A 2 

j£ = W)M» 1 (23) 
in agreement with Fedorov et al., 1992, Eqn (95). For cosmic plasmas, as a rule, A 3> Ro- The helical parameter 
8 — a 2 1 (u 2 ) is generally unknown and hard to reliably estimate; it seems to be rather small for most of the astrophysical 
objects. Below we make the case that the turbulence helicity can, however, be firmly calculated when it is driven by 
twist of the mean magnetic field in the acceleration region, which is likely relevant to particle acceleration in solar flares. 



3 KINETIC HELICITY OF MHD TURBULENCE IN NONPOTENTIAL FLARING LOOPS 

The charged particle acceleration in solar flares occurs in a relatively tenuous corona, where no direct measurement or 
estimate of the kinetic helicity is available. To get a qualitative idea ofpossible origin of the helical turbulence in the 
solar flaring loops and the order of magnitude of the helicity parameter a we consider an MHD turbulence composed 
of random ensemble of Alfven waves with a broad distribution over the wavelength, and, perhaps, containing other 
linear eigen-modes of the magnetized plasma. A single Alfven wave represents in fact an oscillation of the corresponding 
magnetic filed line (or a magnetic flux tube); in the presence of the wave ensemble this field line will experience many 
independent oscillations linearly superimposed on each other until their amplitudes remain small. If the mean magnetic 
field is uniform then no preferable twist is expected in the turbulence composed of these linear wave-modes. If, however, 
the original field line is itself twisted then its oscillations will inherit this twist, resulting in a nonzero helicity in the general 
case of the MHD turbulence generated on top of a nonpotential twisted mean magnetic field, which is, in particular, the 
case of the solar coronal flaring loops. 

To obtain a quantitative estimate we express the helicity parameter in the form a = — r c (/ifc)/3, where hu is the 
kinetic helicity density defined as = u(r,i) ■ V x u(r,t), sf Eqn l[2)l. The advantage of this representation, as we are 
going to demonstrate below, is that in the case of the solar corona we can reliably estimate the kinetic helicity density 
using observational data on the current helicity density (which can be better quantified for the solar corona than the 
magnetic helicity used by Fedorov et al. 1992 although one can be expressed via the other) defined as 

5 A formal quantitative derivation based on the MHD equation set is given in Appendix I. 
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h c = B(r,t)-V xB{r,t) = a FFF B 2 , (24) 

where the second equality is written for a force- free field (FFF), which implies a direct link between the force-free 
parameter lvfff and the current helicity density h c for the case when the magnetic field satisfies the force- free conditions. 
Although none of these two parameters can be directly measured in the corona, we make use of qfff conservation in 
the force-free field along any given field line (B ■ Vqfff = 0) and use the photospheric current helicity measurements 



(or a corresponding nonlinear FFF extrapolation, see, e.g.. iJing et al.l (|2012T ) and references therein), which implies that 



we can have a reliable estimate of the current helicity in the corona. 

Let us now estimate how the turbulence kinetic helicity is related to the current helicity. For the linear MHD modes, 
the fluid velocity is approximately proportional to the magnetic field in the waves, 

u w Tv A b/B, (25) 

thus, assuming p = const for simplicity, we obtain 

v 2 

h k = u(r, t) ■ V X u(r, t) = ^b(r, t) ■ V x b{r, t). (26) 

The last needed step is to estimate the current helicity component related to only the turbulent magnetic field, h c — 
b(r,t) ■ V x b(r,t), which can be done based on the full current helicity h c , Eqn (|24|l . 

To obtain this link let us write down Eqn (|24l) for a field line distorted by an Alfven wave given that even in the 
presence of the wave the force free condition remains on average valid and take into consideration the o?fff conservation 
along this new (distorted) field line, 

<(Bo(r, t) + b(r, i)) • V x (B (r, t) + b(r, t))) = 

a F FF({B (r,t) + b(r,t)) 2 ). (27) 

After averaging Eqn ([3]) , all linear over b terms drop out and then, subtracting the original equality Bo-VxBo = cxfffBq 
from Eqn ([3]), we obtain 

(he) = (b(r,t) ■ V x b(r,t)) = olfff (b 2 ) , (28) 
so the kinetic helicity density takes the form (sf Eqn (|6UI in Appendix I) : 



where (b 2 ) is related to the entire turbulent magnetic field because all waves composing the turbulence are statistically 
independent so their contributions add up incoherently. Now, estimating r c as L c /va, where L c is the main scale of the 
turbulence, see Eqn ([9]), we obtain, in agreement with Eqn (|68[1 of Appendix II, the required kinetic helicity parameter 

Tc(h k ) (b 2 ) 

«~ — % FFF cVA m*' ( ) 

This estimate of the helicity parameter enables us of estimating relative efficiency of the 'helical' acceleration and 
the standard stochastic (Fermi) acceleration given by Eqn (|23|l . which, with the account of Eqns (|30l) and (|25|) . yields 



D h a 2 F FF L 2 c b 2 A 2 2 r2 B 2 (L 



2(*/-l) 

-z r 2 I " tJC \ —2(^ — 1) /oi\ 

D F ^ B 2 Rl - aFFFL ^[lt ) XP ' (31) 

where in the second equality we use Eqn (|10[) to estimate the particle mfp formed by the same turbulence. It is transparent 
from here that the helical acceleration is relatively more important for (i) more nonpotential loops (larger offf), (ii) 
weaker turbulence or stronger mean magnetic field (smaller b 2 /B 2 ratio), and (iii) lower energy (smaller Larmor radius 
Ro). 

We now in the position to make estimates for various nume ric inputs of the relevant pa r ameters. For example , 
substitutinginto Eqn (|31[) modest values of offf ~ 10 -10 cm -1 ( Abramenko et al. 19961 . 1997 ; Longcope et al. 19981 ). 



Ro ~ 1 cm|j L c ~ 10 7 cm, v = 1.5, and (b/Bo) 2 ~ 10 -4 , we find D^/Df ~ 10 s . Note that in flare-productive active 
regions the offf parameter c an be much larger than adopted above: from the nonlinear force free field modeling of active 
region 11158 (February 2011) |jing et~ZI (|2012h reported the mean and maximum values to be |qfff| = 5.2 ■ 10 cm 



and |q!fff| = 5.8 ■ 10 7 cm , respectively. Consider now an extreme case of a strong turbulence (b/Bo) 2 ~ 1, which, 
presumably, happens in a highly twisted flaring loop. T he twist of a quasistationa ry loop, cannot, however, be arbitrarily 



large sinc e the loop is known to become unstable (e.g., Browning fc Priest 19831 ) if the number of twists, Nt w , is large, 



N tw Z 3 (|Aschwandenll2005l . p. 269-270). For a loop with the length L ~ 10 9 cm this translates to qfff ^ olfff.ct 



6 For example, for B = 100 G and T = 1 MK, the thermal gyroradius of electron is ~ 0.2 cm, while of proton is ~ 10 cm. 
© 2012 RAS, MNRAS Q00.nifl6l 



8 Gregory D. Fleishman and Igor N. Toptygin 



4jrN tw /L « 3 • 10~ 8 cm -1 (|Aschwandenll2Q05l . p. 215). For a FFF ~ a FFF , cr we find D h /D F > 1 for L c > 10 5 cm, which 
is likely to be fulfilled for any realistic turbulence model, i.e., acceleration by helical turbulence can be up to a few orders 
of magnitude mo re efficient for low -energy particles in the solar corona than the standard stochastic acceleration by the 
same turbulence. I Jing et al. (2012) reported even one order of magnitude higher values |offf| = 5.8-10 cm implying 
even stronger acceleration, which, however, from the loop stability perspective can only happen in rather short loops, 
L 10 s cm. Note, that in addition to the overall limitation of the helical acceleration energy gain by the value ~ 3eZEhL, 
which cannot exceed 100 keV for any reasonable input parameters, the particle Larmor radius grows with energy, so 
the relative efficiency of the acceleration by helical turbulence drops and the standard stochastic (Fermi) acceleration 
can start to dominate even earlier at some energy < SeZE^L in case of a reasonably strong MHD turbulence. Thus, all 
known advantages of the stochastic acceleration at high energies remain unaffected by the turbulence helicity. 

However, the helical acceleration is efficient at low energies and for a weaker turbulence and so can, perhaps, 
more efficiently, than the standard stochastic acceleration process, extract charged particles from the thermal pool and, 
thus, form a particle seed population from which the particles are then picked up by the main, presumably stochastic, 
acceleration process. Indeed, the turbulence helicity builds up a large-scale regular electric field which can form a runaway 
particle population. T he runaway electrons represent collimated beam s, which are supposed to reveal themselves via 



particle population, i ne runaway electrons represent collimated beam s, wmcn are supposed to reveal tliemselves via 
radio type III bursts ( Aschwanden et alj 19901; Iviahos fc RaoultJ 1995 ) and via sub-second fluctuations of microwave 



( Altvntsev et al hood ) or X-ray (jKiplinger et al.lll983l ; lAschwanden et al.lll996r ) emissions. In particular, if the turbulence 
consists of many turbulent 'cells', where the mean magnetic field has different directions relative to the line of sight, 
then the directions of the electron beams will also be different in agreement with observations in the microwave range 
( Altvntsev et al.l 120071 ). However, at the thermal energy range the Coulomb collisions are important, whose effect we 
investigate in the next section. 



4 IMPLICATIONS FOR PARTICLE SEED POPULATION 

In the considered model of acceleration by helical turbulence, three main ingredients are responsible for the particle 
distribution function formation — Coulomb collisions, stochastic turbulent electric and magnetic fields, and the mean DC 
electric field. Clear division of the particles onto the thermal and nonthermal components in the acceleration region is 
straightforward when binary collisions is the dominant process for the thermal particles (setting the thermal Maxwellian 
distribution), while less important or entirely negligible for high-energy particles (allowing a nonthermal, non-Maxwellian 
distribution). Therefore, to transfer a particle from the thermal pool to the nonthermal component requires a process 
(commonly called 'injection') capable of overcoming the Coulomb losses. We have already established that the role of the 
helical turbulence-induced mean DC field is often more important compared with the stochastic fields at low energies; 
thus, the injected particle population is supposed to be formed in a tradeoff between the (accelerating) DC electric field 
and (decelerating) Coulomb collisions, where the latter can be quantified by the Dreicer field. Note that in addition to 
the Coulomb collisions, the collisionless angular scattering by the turbulent waves can be present even for the thermal and 
suprathermal particles. However, since this scattering does not imply any energy loss, it is inessential to the considered 
here injection threshold and efficiency and so will be neglected, even though it can affect the exact dynamics of this 
process. To be specific, we adopt that there are no turbulent waves capable of resonantly interact with the thermal and 
slightly suprathermal particles because the corresponding linear modes have a strong damping; thus, the thermal particle 
transport is mediated by the Coulomb collisions. Let us estimate how big the regular electric field induced by the helical 
turbulence can be compared with the electron Dreicer field. Substituting Eqn (|30[) into Eqn (|15|l we obtain the electric 
field in the form 



/ 2.4 -10- 7 \(_B_V( n e 
h y 7.3 • KT 5 J \W0 GJ Vl0 10 cm-3/ 

/ a FFF \ ( L c \ f(b 2 )/B 2 \ ( statvolt/cm \ 
V10- 10 cm- 1 / V10 8 cm) \ lO" 4 ) I V/cm J' 

which, for reasonably modest assumptions on the acceleration region parameters yields 



E h ~ 10~ 5 - 10~ 4 V/cm. (33) 

Efficiency of particle acceleration by a DC electric field is specified by its ratio to the Dreicer field, whose value is defined 
by the particle Coulomb collisions. In particular, the electron Dreicer field is 
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elnAc I 2-l(T 7 

Eoe = o " 



6- 10 _J / VlO 10 cm- 3 



T \ /In Ac \ ( statvolt /cm 



10 7 KJ V 20 J \ V/cm 



(34) 

therefore, in our example the induced electric field, Eqn (|33|l . is slightly sub-Dreicer, although it can exceed the Dreicer 
field, e.g., for more tenuous plasma and/or for stronger turbulence. The sub-Dreicer field creates beams of the runaway 
electrons with the velocity exceeding the critical velocity, v c ~ VTsy/Eoe/ E, which can immediately produce the type 
III radio bursts and then be picked up by the main acceleration process, i.e., the seed electron population is formed 
by the runaway electrons, whose amount and, accordingly, the thermal-to-nonthermal energy partition, depends on 
the E/Ede ratio. In case of super-Dreicer field, E/Edc Si 1, most of available electrons will 'run away' and form the 
seed population, i.e., almo st all these electron s are eventually acceler ated in agreement with recent observations of the 
flare acceleration regions (jrCrucker et al.lboid : iFleishman et alj|201ll ). Thus, we see that the helical component of the 



turbulence is easily capable of forming a seed population of electrons, which is needed for the mechanism of stochastic 
acceleration to efficiently work. 

Consider now how this electric field can affect formation of the ion seed populations. For the estimate let us adopt a 
three-component plasma consisting of electrons (e), protons (p), and one more sort of ions (i) with the charge Z\e\, mass 
rrii, and number density m so that n e = n v + Zrii, permitted by a sub-Dreicer uniform electric field E. To address the 
question of the seed population formation we have to estimate the value and direction of the component flow velocities in 
the given plasma. To do so we have to determine conditions providing the balance between the electric force and forces 
of the dynamic friction between the plasma components driven by the Coulomb collisions between the particles. 

The dynamic friction force acting from a plasma component a on a given 'te st' particle is ca lculated from the 
averaged momentum exchange between this plasma component and the test particle l|Trubnikovll 19651 ) . i.e., 



T a {u) = f -^—!L-f a ( v) d\ (35) 



where u is the velocity of the test particle, ii a = Mm a /(M + m a ) is the reduced mass defined by the test particle mass M 
and a-component particle mass m a , f a (v) is the distribution function of the component a, and Q = 47re 4 lnAc. To find 
the force acting on a 'mean' test particle we have yet to average force (|35[) over distribution function of the 'test' particle 
component. Both these averagings are convenient to perform assuming u — ut + Su and v = vt + Sv, where uj and vt 
are the corresponding thermal velocities with zero means (but nonzero rms values), while Su and Sv are flow velocities 
driven by external electric field. Since we adopted E/Eoe < 1, the denominators are defined by the thermal velocities, 
whose contributions to the numerators are zeros, so the numerators are solely determined by the flow velocities. 

Accordingly, the balance of forces acting on the electron component yields (we drop 5 and use v for all Sv below for 
short) 

-\ e \E-J£±-(v a -v p )-^£-{v*-Vi) = (36) 
where we took into account that vts^>vt p , VTe^vn and fx a ~Tn e . Then, in a similar way, we write for protons 

le\E--Q?l-(v p -v e )--?^(v p -v i )=0, (37) 

where we adopted for simplicity VTp^vn; this approximation is the least accurate for the Helium ions but even is this 
case the error is within 30%, which is acceptable for the purpose of the estimate; with the same accuracy we use below 
the proton mass m p for the reduced ion-proton masses rrn p . Then, we add up Eqns (|36[) and (|37[) to eliminate the electric 
field: 

QZ{Z - l)m Z 2 Qn t Z 2 Q ni 

=— * — v e =— v p H =— Vi = 0, (38) 

rn e ?4 e m tp v^ p m lp vj, p 

where the plasma quasineutrality n e = n p + Zrii is taken into account and some small terms are discarded, and then 
eliminate the electron velocity using the momentum conservation law v e = — (n p m p Vp + riimiVi) / (n e m e ) (provided that 
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Figure 1. Illustration to the origin of the heavy ion velocity directed oppositely to the external electric field due to dynamic friction 
from moving electrons. 



the plasma is immobile as a whole in this reference frame Jj 



1 + 



mi p v Tp Z — 1 rumi 



m e n e 



m ip v Tp Z — 1 n p m p 



(39) 



This equality demonstrates a remarkable property of the admixture ion behavior in the external electric field: the 
admixture ions move in the direction opposite to the main ions (since the second, negative term dominates the rhs for 
the natural abundance of the protons), i.e., the admixture ions move to gether with elect r ons in the direction opposit e 
to the electric vector direction. This anomalous behavior noted earlier bv lGurevichl (|1961T ); iFurth fc Rutherford! (119721 ) ; 
Holman 1 19951 ) originates from the fact that the dynamic friction force produced by moving electron component on the 
ions with Z > 1, and so proportional to Z 2 , turns out to overcome the electric force, which is proportional to Z. Another 
remarkable property revealed by Eqn (|39[) is that a significant ion flow velocity can only be achieved for a relatively 
tiny ion population for which the second term in the lhs brackets can be neglected compared with 1; for some abundant 
ions including helium and oxygen it is not the case since the second term is larger than or comparable to 1 even though 
ni/n p <C 1 for them. 



Let us make estimates of the flow velocity for the Helium isotopes. Consider 4 He [Z — 2 and rm 
and adopt mv% ~ T for all plasma components, i.e., VTp/vTe ~ y/ m e /m p , then 



1 + 2- 



2n e 



7Ti4 = 4m p ) first 



(40) 



Since ri4 /n e ~ 0.1 > V"W (4m p ) ~ 1/86 ~ 0.012, for a rough estimate we can discard the first terms (ones) in 
both brackets, which simply yields V4 ~ — (n p /Arn)v v . To express the ion flow velocities via the electron flow velocity, 
however, we have to calculate the ion velocity more accurately, i.e., retain small terms in Eqn (|40[) . which yields: 



v 4 « - 



4n 4 



1 - 



n e (4ri4 + n p ) lm e 
2n p ri4 



(41) 



' We note that in a purely electron-proton plasma the momentum conservation would imply a very small proton flow velocity 
v p = —(m e /m p )v e . 
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Substituting Eqn (|41[) into the component momentum conservation law we find for the protons 

2n 4 [ml . . 

v p ~ n — \/ — Ve ' ( 42 ) 

n p + 4ri4 y m p 



which is valid if 2n4,/(n p +4714) > \J m e /m p and, accordingly, it is noticeably larger than v p — — (m e /m p )v e in the 
purely hydrogen plasma. Now, substituting Eqn (|42p into Eqn (|41|) . we find for the 4 He ions 

Tip / TTl e . . 

f 4 ~ —7 £ — ; r-,/ V e . (43) 

2(n p + 4n 4 ) y m p 

A similar analysis for the other helium isotope, He (Z = 2 and uij = m3 = 3m p ) gives rise to essentially different 
result because of much lower 3 He number density, 713 ~ 1.7 • 10 _4 ri4. For such a low density, we can safely discard the 
second term in the lhs of Eqn (|39|) (while still discard 1 compared with the second term in the rhs), which eventually 
yields: 

V 3 « , , . rV e , (44) 

n e (n p + 4n 4 ) 

roughly one order of magnitude larger than the 4 He flow velocity V4, Eqn (|43[1 . 

The results obtained allow estimating 'effective Dreicer' fields for various ions, which we define as an electric field 
needed to form the ion flow velocity equal to the thermal velocity of the same ion. Rewriting Eqns (|43p and (|44|l in 



a compact form v% = KiV e and taking into account that v e — —VTeE/ Edb and VTe = vn y rru /m e (we still assume 
equal temperatures of all plasma components), we obtain Vi/vn = —Ki^/mi/m e E/EDe = —E/Ed%, where the latter 
equivalence is a definition of the effective ion Dreicer field. Then, using Eqn (|42[) for protons and Eqns (|43|) and (|44|) for 
Helium ions we find 

E Dp ^^±l±E De , (45) 
gpaM M!V±±u) .& g ., (46) 



n p ri4 y 7713 

Sd4 „ ^±i^i £De . (47) 

Tip 

Adopting ri4 = 0.1n p (i.e., n e = 1.2n p ), we find Eds ~ 0.23i5o e , while _E_d4 ~ lAEoe and ~ 6i5c e , i.e., Edz^Edc < 

EdA <C £-Dp- 

This remarkable difference in the effective Dreicer fields implies that a given electric field can be sub-Dreicer for the 
electrons and ions with large abundance (e.g., p and 4 He), while super-Dreicer for ions with small abundance (e.g., 3 He), 
which is in particular the case for the electric field driven by the helical turbulence estimated by Eqn (|33|) . Therefore, 
all available 3 He ions can 'run away' in the regular electric field present in the helical turbulence and, thus, all the 3 He 
ions become available for further acceleration by this turbulence, while only a minor fraction of the 4 He ions from the 
corresponding Maxwellian tail will run away in the same electric field, so a much smaller number of the 4 He ions are 
available for acceleration, which offers an efficient mechanism of 3 He enrichment in solar flares. Note, that given that 
Ed4 <^ Ep v , the same mecha nism can be responsible for the 4 He/p enhancements discovered from the gamma-ray line 
study jShare fc Murphvlll99g| ). 



5 DISCUSSION 

Currently, it is widely accepted that the stochastic mechanism plays a key role in charged particle acceleration in solar 
flares. Any stochastic acceleration model implies a turbulent spectrum of waves (linear or nonlinear eigen-modes), whose 
resonant or nonresonant interactions with the charged particles result in the particle energy gain on average. In this study 
we have found that the linear MHD modes excited on top of a twisted, nonpotential mean magnetic field necessarily 
possess a nonzero kinetic helicity and so differ from the standard MHD modes excited on top of the uniform mean 
magnetic field. Consequently, the turbulence composed of these 'modified' MHD modes will also possess a nonzero 
kinetic helicity. 

We, therefore, have established that nonpotentiality of the magnetic field at the flare site (i.e., twists of the magnetic 
field lines) implies necessarily a certain level of helicity of the MHD turbulence created there by the primary flaring 
energy release. This turbulence helicity has potentially a number of highly important effects on the particle acceleration 
in the flares because, in addition to the common for the turbulence random magnetic and electric fields, capable of 
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particle scattering and stochastic acceleration respectively, it also generates a nonzero mean regular DC electric field 
with a rather significant value, see Eqn Q. 

For a plasma obeying the standard Ohm's law with Spitzer conductivity such an electric field would represent 
a problem since this would imply an unrealistically strong electric current, J c ~ 10 18 A, and correspondingly strong 



magnetic field, B c ~ 10 G (e.g., lEmslie fc Henouxl 1 19951 ). For the considered here helical turbulence the situation 



is entirely different. Indeed, the enhanced turbulence reduces the plasma conductivity (by introducing an anomalous 
resistivity) and modifies the Ohm's law accordingly, see Appendix II, and so the role of the enhanced, helicity-driven 
DC electric field is just to support the original electric current needed to maintain the original magnetic field twist 
(jocVx B), while not to enhance the current or generate additional magnetic field. 

This electric current is formed by steadily moving thermal particles for which the effective electric conductivity 
derived in Appendix II applies. However, because the electric field is enhanced, the amount of runaway electrons capable 
of leaving the thermal pool and transiting to a nonthermal population increases accordingly. We emphasize that the 
helicity-driven DC electric field has almost no effect on particles with energy higher than 10s keV, where standard 
stochastic acceleration remains dominant. Accordingly, the particle spectra at <^ 10 keV will presumably have a standard 
power-law form as expected for the standard stochastic acceleration process. In contrast, the DC field can affect and even 
dominate the particle acceleration at low energies — from the Maxwellian tail to a few or a few 10s keV, where the rate of 
stochastic particle acceleration can be greatly enhanced when the turbulence is helical, i.e. when the helicity parameter 
a^O; see estimates of Eqn (|31|) . 

There, depending on how strong the DC electric field is relative to the Dreicer field, it can control what fraction of the 
plasma particles will 'runaway', forming a seed particle population, and so go eventually into the accelerated component 
due to the standard stochastic acceleration of the seed particles. This process offers an efficient mechanism of the flaring 
energy partition betw een the thermal plasma and accelerated nonthermal particles demanded by observations (see, e.g. 



Altvntsev et al.l 2012 ). which is in principle capable of providing any partition — from 'heating without acceleration' 



( Battaglia et al.l 120091 ') to 'acceleration without heating' (jFleishman et al.l 120111 ) . These runaway electron populations 



can form beams directed along the mean DC electric field, capable, in some cases, of producing type III radio bursts 
and some of rapid fluctuations observed from the microwave and HXR radiations. In fact, under conditions of angular 
scattering a beam can only exist over roughly one mfp, since the beam-forming electrons will then isotropize. Therefore, 
if the beam produces or does not produce a type III burst depends on how the mfp A relates to the density scale, Ld'- 
a drifting type III burst passing a broad range of frequencies can only be produced if A S> Ld, while the beam will 
be destroyed before a significant spectral drift is achieved if A <C Ld, thus, the spectral range of the corresponding 
radio emission must be rather narrow in the latter case, leading to a 'f ailed type III burst ' . Perhaps, such fa iled type III 
bursts represent a widely observed phenomenon called 'metric spikes[f| ijBastian et al.lfl~998l : llBenz et al.lll996l ). which were 
proposed to come from the acceleration region directly. Indeed, if the level of turbulence is enhanced in the acceleration 
region and the mfp is reduced accordingly, the beams will be destroyed and the beam-driven instability responsible for 
the coherent radio emission quenched shortly after the beam formation; thus, the elementary radio burst will remain 
narrowband in contrast to formation of the broadband drifting structure in the classical type III bursts, which originate 
outside the acceleration region where the mfp is noticeably larger. 

Furthermore, we point out that this same electric field can form highly non-even seed populations of various ions 
because of correspondingly distinct effective Dreicer fields of the ions with different abundances, which offers a natural 
way of understanding of the accelerated particle population enrichment by ions with relatively small abundances. This 
enrichment mechanism based on electrically selective extraction of the ions from their original thermal distributions 
with natural abundances, described by Eqn (|39[) . can account for t he entire 3 He e nhancement in the 3 i/e-rich events, 
which have been puzzling for about 40 years by now (see, however, Liu et al. 200rJ ). Eqn (|39[) shows that the ions with 



Z > 1 and a low number density can have relatively small effective Dreicer fields, which is favorable for them to run 
away and be picked up by the bulk acceleration process. In contrast, no enrichment is predicted for the 2 H ions, which 
are a comparably rare as the 3 He ions, but have Z — 1 and so, for the hydrogen isotopes 2 H and 3 H, the terms 



containing Z — 1 drop out from Eqn (|39[) and these ions have exactly the same flow velocity as the protons implying no 
enrich ment for any of the hydrogen isotopes. And indeed, no 2 H or 3 H enrichment has ever been observed from solar 
flares ( Kocharov fc Kocharov 19841 ) in agreement with th e proposed her e non-resonant enrichment mechanism. Note, 



that although resonant enrichment mechanisms (see, e.g.. lLiu et al .1120061 ) are also capable of producing He and other 



isotope enhancements, it is yet unclear how they are always so fine tuned to keep the 2 H ions unenriched in all cases. 

On the other hand, high-Z ions can display even stronger enrichment than the 3 He ions, because (Z — 1)/Z ~ 1 
for them in contrast to helium, where [Z — 1)/Z = 1/2. In fact, enrichment of 3 _ffe-rich events with ultra-heavy ions is 



8 The metric spikes are narrowband spikes originati ng at 200-400 MHz at the frequency range adjacent to the starting frequency 
of the classical type III burst (see, e.g., Figure 10 in lBastian et al1ll998l) . 
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widely observ ed 1 Masonl 120071 ). Note also that a high correlation between the 3 _ffe-rich events and radio type III bursts 
I Mason 2007) also receives a natural interpretation within the proposed enrichment mechanism, because the formation 
of the 3 He seed population happens along with formation of the runaway electrons responsible for generation of the radio 
type III bursts. As has been said, whether the electron beam is formed or not is determined by the balance between the 
runaway electron acceleration by the electric field, which is proportional to A a B^/(b 2 } and their angular scattering by 
the same turbulence, which is reciprocal to A. Therefore, for a stronger turbulence, when, e.g., the induced electric field 
becomes closer to the 4 He Dreicer field, Eqn (|47[) . the electron angular scattering can effectively isotropize the electron 
beams and so quench the type III burst generation in more powerful flares not showing a considerable 3 He enrichment 
I Mason 2007), although type III bursts can still be produced outside the acceleration region due to electron transport 
effects. 

Then, in the presence of this electric field the particle flux direction, according to Eqn (|19b[) . depends now on their 
electric charge sign; thus, the accelerated electrons and protons are supposed to precipitate into the opposite foots of 
the same flaring loop or loop system, which offers a consistent explanation of the observe d spatial displacement between 



sources of HXR (produc ed by electrons) and gamma-ray (produced by ions) emissions jLin et alj [2003 ; iHurford et al 
200rj ; IVilmer et al.ll201lh 



6 CONCLUSIONS 

We have argued that the presence of strong current helicity (twists of magnetic filed lines) implies a proportionally strong 
helicity of the turbulence produced by release of the flaring energy stored in the helical (nonpotential) magnetic field, 
which, in its turn, has highly important implications for the stochastic particle acceleration in solar flares. The proposed 
mechanism offers a nice consistent way of interpretation for (i) thermal to nonthermal energy partition in flares, (ii) 
electron beam formation, (iii) preferential acceleration of low-abundance isotopes with Z > 1, and (iv) observed spatial 
separation of the electron and proton precipitation sites. Therefore, the stochastic acceleration by helical turbulence is a 
promising mechanism in the magnetized solar corona including the solar flares, jets, and eruptions. 

It should be noted that this attractive mechanism of particle acceleration by helical turbulence has not yet been 
analyzed in any detail, so many important aspects of it remain unclear. Realistic self-consistent nonlinear models of 
particle acceleration by the helical turbulence supplied by (yet nonexistent) nonlinear theory of the helical turbulence 
are called for to fully assess the role of this acceleration mechanism. 
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APPENDIX I 

HELICITY OF MHD EIGEN-MODES IN A PLASMA PERMITTED BY A NONPOTENTIAL 
MAGNETIC FIELD 

Let us consider the kinetic helicity of the linear MHD modes and corresponding mean electric field Eh — — (it x b) /c in 
a plasma permitted by a nonpotential (and constant in time) force-free mean magnetic field based on the original MHD 
equations for the magnetic field and fluid velocity: 

V B = 0, 

dB 

-— = V x [u x B] + u m AB, v m = - — = const, (48) 
at Ana 

p + (u-V)uJ = -Vp + f+^[VxB]xB + T t Au + |v(V-it), (49) 

where v m is the magnetic diffusivity, a is the electric conductivity, p and p axe the plasma density and kinetic pressure, 
/ is an external force (e.g., gravity), and rj is the dynamic viscosity. To investigate properties of weakly damping linear 
modes we neglect the dissipative terms entirely. Then, for simplicity we only consider the helicity originating from the 
magnetic field twist (the corresponding pseudoscalar is formed by the dot product of the polar vector j and the axial 
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vector B: j ■ B, while no other polar vectors are explicitly considered) and so we neglect the kinetic pressure and the 
external force assuming that their contribution to the helicity is smaller; so the equations read: 

VB = 0, ^=Vx[uxB], (50) 

^ + (u-V)u) = ^[VxB]xB. (51) 

Now, to describe the small-amplitude linear modes satisfying these equations, we have to linearize them using 
Eqn (|14|l for the magnetic field and adopting that b and u to be the first-order oscillating values of an MHD mode. 
Importantly, that upon substitution of Eqn (|14|l into Eqn (|51[) we have to take into account that in the nonpotential 
force- free field V x Bo = ufffBq ^ unlike cases of a uniform field or nonuniform potential field, which yields 

§ = (Bo-V)u-(u-V)Bo, § = ^BoxH-^(Vxb)xB». (52) 
at at vkq inp 

Within the eikonal approximation (i.e., adopting the wavelengths of the eigen- modes to be small compared with 
the source inhomogeneity scale) we can write b = fo^e 1 *^ -1 "' and a similar for u, which yields equations for the 
corresponding complex amplitudes: 

b = (i/u)[i(B ■ Vip)u - (u ■ V)B ], (53) 
u = ^^Bo xb- -r^—[(b • B )V^ - (Bo ■ VV>)6]. (54) 

inpuj i%pU) 

Let us solve these equations to the first order accuracy over the small parameter olfff /|V^| < 1. In the zeroth 
order approximation we have 

b=--(B -Vi^)u, u = --^—(B - Vi>)b, (55) 
ui inpoj 

which yields the eikonal 

V\\ip = ±oj/va, va = Bo/^ATip. (56) 

It is easy to see that in the zeroth approximation these perturbations are identical to the purely alfvenic modes for which 
the conditions b ■ Vi/> = 0, u ■ = 0, b ■ Bo = 0, and u x b = are satisfied. Since we use the complex amplitudes the 
bilinear terms like ab must be computed as (l/2)K(afo*), where in addition to averaging over the period T — 2n/ui we 
also average over the random phases of Fourier harmonics: 

M;> = (i/2xO(tf/«,-44). (57) 

In the zeroth over cifff approximation, the kinetic helicity parameter is zero: {u • (V X u*)} = 0. 
Now, taking into account the first-order over \o.fff /V^>| terms, Eqn (|54|) yields 

(Bo ■ VVOb + — — [iuFFFBo x b + (6 ■ B o )V<0], (58) 



AllpUJ 4-7T pUJ 

where b ■ Bo 7^ and 

V x u = -— — (Bo ■ V-0)(VV x b) + -7— — [ — Vt/ 1 x (b ■ V)B + b x (B ■ V)W + b x (Vib ■ V)B 1. (59) 
Anpuj Aivpui 

These two expressions allow calculating the kinetic helicity density 

isft(t t -(Vxn*)> = ia FFF ^^ ) 

Z 4 S3, 



{h k ) = -$l(u-(Vxu*)) = 1 aFFFV A ±£f, (60) 



which is nonzero any longer; it coincides with estimate (|29[) in the order of magnitude. 



APPENDIX II 

MEAN DC ELECTRIC FIELD INDUCED BY THE HELICAL TURBULENCE 

Let us now calculate the large-scale electric field, created by helical turbulence, taking into account the entire range of 
harmonics forming the random fields of u(r,t) and b(r, t). From the first of Eqns (|52l) we find 

b{r,t) = Bo, f du ( r ^) dT _9Bo j* u (p ) d (61) 
J-00 ox v dxp J_ oa 

which allows to form the required bilinear cross product u x b and perform the averaging 

E h = -(l/c){u{r,t)xb(r,t)) (62) 

© 2012 RAS, MNRAS 000. MTE\ 
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over the turbulent ensemble according to Eqn (|15l) . which is convenient to express via the vector components 

(u a (r,t)b ll (r,t)) = -^ f* {u a (r,t)u (r,T)}dr + B Ov f (u a (r,t)^^)dr. (63) 

The first term in the rhs describes a modification to the Ohm's law due to the wave turbulence ensemble (anomalous 
resistivity), while the second one describes another correction to the Ohm's law due to the turbulence helicity; we will 
see below that these two terms provide comparable contribution in the problem studied. Adopting for simplicity that the 
turbulence is statistically uniform and isotropic, we can express the first integral using the invariant Kronecker's tensor 
as 



/' 

j — e 



(u t7 (r,t)uf i (r,T))dT = ((« 2 )r /3)<5 ai 9, (64) 



where (u 2 )r c /S — v t is the magnetic turbulent diffusivity. 

The second integral, which is a third rank tensor, is obviously proportional to the Levi-Civita's permutation tensor 
e CTM i, and the kinetic helicity pseudoscalar a 



8u u (r,r)\ dT= l 



u a (r,t) — — )dr=-ae, fw . (65) 

Now, the large-scale DC electric field can be written in the form 

E h = (i/ t /c)V x Bo - (a/c)B = -(a eff /c)B , (66) 
where, according to Eqns ([60]l and ((2| 



and 



12 -,„.- A ir, a R ff = a-{vt/c)a F FF (67) 



is the total effective helicity parameter specifying the electric field in a nonpotential flux tube. Although both terms 
in a e ff are comparable, the contribution from the anomalous magnetic diffusivity (anomalous resistivity) is four times 
larger than the direct helicity contribution. The correlation time can be estimated using the correlation scale L c and 
the Alfven velocity va- t c = L c /va- Finally, recalling that the magnetic field vector can be expressed via the electric 
current density, Eqn (|66[) can be presented in the Ohm's law form j = a e ffE, where <J e ff = —(c 2 /4tt)(ocfff /a e ff) is the 
anomalous (turbulent) electric conductivity. When the effective conductivity is much smaller than the collisional Spitzer 
conductivity, having the same electric current (i.e., the same magnetic twist) requires a proportionally larger electric 
field, which is induced by the turbulence helicity. 
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